Thermodynamic origin of universal fluctuations and two-power laws 
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We discuss universality of response functions in systems 
with excited degrees of freedom. We propose a unification of 
two existing phenomenologies, two-power law decay and de- 
viation from power law due to non-extensivity. A universal 
curve is derived by maximizing entropy with a non-linear con- 
straint. The same formalism can explain the universal fluctua- 
tion curves which have been discovered recently by Bramwell, 
Holdsworth, and Pinton. 

PACS numbers: 31.70.Hq, 05.40.-a, 05.70.Ce 



Many experiments measure the probability p{t)dt that 
an event takes place within the time interval [t, t -\- dt] 
after excitation of the system at t = 0. Classical argu- 
ments predict exponential decay. We are interested here 
in a power law decay at large times 
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with < a < 1. 

There is an extensive literature explaining non- 
exponential decay by means of fractal theory. In this 
domain of research it is accepted that a scaling law like 
(|l|) holds only for a limited range ti < t < t2 oi the argu- 
ment t. For t < ti collective effects presumably restore 
non-fractal behavior. For large t individual behavior on 
microscopic scale dominates and non-universal decay is 
expected. A nearly ideal example of this behavior has 
been observed recently in quantum dots |^ , with a power 
law extending over a range of 10^ of time scales. Many 
theoretical models explain power law decay as the correct 
asymptotics in the large time limit. 

There is more and more evidence that also non-power 
law decay can be described by a universal probability 
density function (PDF). An early explanation Q of non- 
exponential decay involves the notion of an activation 
energy required to trigger the decay process. These au- 
thors proposed the following modification of (|^) (using 
our notations) 



p{t) - (l + a(l-a)cc;t)-i/(i-") 
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with a some positive constant, and uj a constant with di- 
mension of a frequency, introduced to obtain the dimen- 
sionless combination cut. Later on, (H) lost importance in 
favor of stretched exponential decay H 



which is expected in case of processes with a scaling dis- 
tribution of relaxation rates. Here, we restrict ourselves 
to the case of asymptotic decay with a single power law. 

In dielectric relaxation one observes usually two re- 
gions of power law response (in our notations) 
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(a)i)-l/(l-") ujt < 1 
(a)t)-l/(l-'') Lut:$>l. 
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(see e.g. Weron and Jurlewicz Q,^). The point of view 
of the present paper is that both behaviors, as described 
by (H) and (|^), should be combined. More precisely, the 
generic behavior in the region cut ^ I should be described 
by (||) instead of (|^). This substitution makes sense be- 
cause (H) describes a curve which starts off linearly and 
bends over towards a power law decay with exponent 
-1/(1 -a). 
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FIG. 1. Log-log-plot of p{t) with a — 500 for different com- 
binations of exponents (a,p): (0.1,0.1) dotted, (0.1,1) solid, 
(0.1,0.8) short-dashed, (0.1,1.08) long-dashed. 

The actual formula, which we will derive below, and 
which covers the whole fractal region, is (see Fig. 1) 



pit) - (l - a + a (1 + (1 - p)u}t) i-p ) 



-1/(1-") 
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Three regions can be distinguished. In regions 1 and 2 the 
decay is as described by (0) . Regions 2 and 3 show power 
law decay according to (0). Note that p > 1 corresponds 
to a super-exponential decay. 

For p = a expression (|^) reduces to (^). If a = 1 then 
it reduces to (^ with a replaced by p. In the limit p = 1 
one obtains 
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p{t) ^ (1 — a + aexp((l — a)ujt)) 



In the limit a = 1 one obtains 



p{t) ^ {1 + (1 - p)ut) 



-a/{l-p) 
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In any case, in the hmit a = p = 1 the decay is exponen- 
tial p{t) ~ exp{—aLut). 

We found first evidence for (^) in the work of Tsallis, 
Bemski, and Mendes They reanalyze the experimen- 
tal data of ||] and observe that the description using (||) 
improves significantly by using formula which was 
introduced on an ad hoc basis. The full data set of 
could be described with it. The universality of (B) is 
further supported by recent work of Montemurro |M] us- 
ing a very convincing analysis of linguistic data. In both 
papers |^ and |^ there is evidence that the asymptotic 
behavior for large t still deviates from and follows a 
power law instead of decaying exponentially. In the logic 
of the present paper, this is a cross-over from one power 
law to another, as described by (Q). Hence (H) should be 
used with p ^ I instead of (^) . In an ad hoc modifica- 
tion of (^) was proposed, slightly different from (^). The 
same expression was used in and fits the experimental 
data over the whole range of t-values. It is clear that 
will fit as well, but has a number of advantages: 1) 
there is a general formalism, which is used to derive it; 
2) the PDF has a workable algebraic expression, which 
is not the case for the formula proposed in |^ ; 3) within 
the context of the formalism an average time {{(^t)) is 
defined using a nonlinear average - see (|l|) below. The 
average time characterizes the PDF p{t) together with 
the two exponents a and p, and the cross-over frequency 
uj. It replaces the fitting parameter a which has a less 
clear physical meaning (see Fig. 2). In particular, in the 
case p = 1, we show that ((•)) satisfies a property of 
additivity, which implies that it is meaningful to add av- 
erage times of different experiments all having the same 
values of a and uj. 
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FIG. 2. Kolmogorov-Nagumo average time {{ijjt}) as a 
function of 1/a for several values of (a,p): (0.9999,0.9999) 
straight line, short-dashed; (0.999, 0.9999) long-dashed; 
(0.999, 1) solid line; (0.8, 0.9) thick-dotted; (0.9, 0.9999) 
thin-dotted. 



Up to now we considered decay as a function of time. 
However the PDF (|) can also be considered as a general- 
ization of the distribution of Boltzmann-Gibbs. Time t is 
then replaced by energy frequency w by a fixed inverse 
temperature /3. With this interpretation (^) is the equi- 
librium distribution of non-extensive thermostatistics . 
The entropic parameter q of this formalism, coincides 
with the exponent a (the relation becomes q — 2 — a in 
the recently modified formalism ||l^). In the context of 
non-extensive thermostatistics ( p|) h as numerous applica- 
tions — for a review see Q or |12{ . In this context (||) 
becomes the equilibrium distribution for thermostatistics 
based on Renyi's choice of Komogorov-Nagumo averages 
and a-entropies - see . 

Let us now discuss universal distributions of fiuctuat- 
ing quantities. Their universality has been discovered 
recently by Bramwell, Holdsworth, and Pinton (BHP) 
[|l3|,0, based on earlier work |p^. The result received 
ample support in the literature |l6|-p^. The BHP den- 
sity function is of the form 



3BHp(e) ^ exp (%- 
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with y = c{e — u). Setting 6=1 one finds the famous 
Fisher-Tippett density appearing in statistics of extremes 
(see Gumbel |^). However, the value which empirically 
best describes fiuctuation spectra is & = tt/2. 

Using the standard central limit theorem (CLT) we 
show at the end of this Letter that our formalism implies 
(for p > 2/3) 
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g{e) = ( expp [v + e] ^ 

^{-^['Papiv + e) - (l)ap{v)]y 
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with CT^ the variance of p{t), g{x) the normal density, 
and the functions (j)ap and exp^ defined later in the text 
(for p < 2/3 the normal density function g(e) should be 
replaced by a Levy-stable PDF). 

Of particular interest is the limit p — I (corresponding 
to Renyi a-entropies - see further on). Then (^ can be 
written as 
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with y = (1 - a)e and d = exp((l - a)v)/N'^/'^a{l - 
a). This expression differs in an essential way from (|^). 
However, both PDFs have a similar shape (see Fig. 3). 
Note that (^) has an additional fitting parameter, the 
role of which is taken over here by the exponent p. Very 
accurate data will be needed to distinguish the two PDFs 
on an experimental basis. An advantage of ( p^ is that 
in the limit a = 1 it reduces to the normal density with 
width (T, as expected. Note also that, if .g(e) were replaced 
by the exponential density, and p = 1, then (0) reduces 
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to (^) with 6=1 and c = I — a. However, the exponential 
density is not Levy-stable. Therefore we do not expect 
it on thermodynamic grounds. 
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FIG. 3. Comparison of g{e) and (7bhp(^)- 9i^) with p — 0.8 
dotted, with p = 0.9 short-dashed, with p — 1 long-dashed; 
Sbhp(^) solid; the curves are normalized to 1 in the origin; 
parameters of g{e) are v — 1, a = 0.5, and = 4N; the 
parameters of g^Hp (e) have been adapted to fit roughly g{e) 
with p = 1. 

In the final part of this Letter we derive formulas ^ 
and using thermodynamic arguments. Weron and 
Jurlewicz express the response function p{t) in terms 
of PDFs Pi{t), describing relaxation of individual dipoles. 
To do this, they use statistics of extremes |25j. The un- 
derlying assumption is that the whole system relaxes as 
soon as one of the individual dipoles relaxes. Our as- 
sumptions are not so drastic. We start with a thermo- 
dynamic argument, and replace statistics of extremes by 
the usual CLT. 

About 20 years ago MontroU and Shlesinger pro- 
posed a thermodynamic derivation of Einstein's diffusion 
law. It follows by optimizing entropy of the PDF p{t) 
under the constraint of zero first moment and fixed sec- 
ond moment. They observed that anomalous diffusion 
can be explained along similar lines using an ad hoc con- 
straint replacing the constraint on the second moment. 
Alemany and Zanette |2^,^, followed by Tsallis et al 
|^,|o|, repeat these arguments, replacing Shannon's en- 
tropy by the generalized entropy used in nonextensive 
thermodynamics H 



1 - a 



(l + (l-p)x)(^-")/(i-'') -ll . (13) 



The a-deformed exponential and logarithmic functions 
are defined by p3\ 



exp„(.T) = [1 + (1 - a)x] 
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The inverse function is ippaix). Following 
sponding definition of entropy is 



(14) 
the corre- 



Sap{p) = (t>p 



dtp{t) ln„(l/p(t)) 



In the limit p = 1 this is the a-entropy of Renyi 



^ai(p) = Y^ln( / dtp{tr 
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On the other hand, one has (l)aa{x) = x so that Saa{p) 
coincides with (|l^). Let us now optimize entropy Sapip) 
under the constraint that the average decay time {{(^t)) 
has a given value. A straightforward calculation using 
Lagrange parameters produces an implicit expression for 
the PDF p{t). By introducing a free parameter a the 
latter can be made explicit, resulting in (||) (see |^). 

In order to calculate the fluctuation spectrum assume 
that p > 2/3 (then the second moment of p{t) exists 
and CLT holds). If A{t) is an arbitrary random variable 
we will denote its linear (experimental and theoretical) 
averages by 



N 

(A)A. = ^^(tfe)/Ar, {A}= / 
fc=i 



dtp{t)A{t) 



Denote by P{x) the probability of an event x, given the 
PDF p{t). In this notation CLT can be formulated as 



P 



[{A)^<z)^F[^iz~{A)) 
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with (T^ the variance of p(i), F the normal distribution 
function, and N sufficiently large 
The nonlinear averages 



So^ip) = / dt 
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a>0,a^l. (11) {{A))^^r\{m)N), {{A))^r\{m)) 



These authors introduce a suitably adapted constraint 
on the second moment. Recently the present authors 
generalized nonextensive thermostatistics by founding it 
on nonlinear Kolmogorov-Nagumo averages l3lU32] 



((/» = <P- 



dtp{t)4>{f{t)) 



(12) 



This average depends on a monotonically increasing func- 

with 



tion 



x). In the present paper we choose 
<ap{x) = lnQ(expp(x)) 



-'ap 



are related via 4>~^ to the linear ones of the random vari- 
able B{t) = 4)[A{t)^ . Since </) is monotonically increasing 
we can write 

p{{{A))n < {{A)) +e)= P(WA))^ < cbirn) + e)) 
and apply CLT to its right side. One finds 

p{{{A))^<{{A))+e) 

«^(^W((^))+^)"<^(((^))))- (18) 
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Introduce a limiting PDF ^(e) by 



-p[{{A))^<{{A))+e)^^9{e). (19) 

(j)ap we finally 



Using g(x) = dF{x)/dx and choosing 
obtain 



K6) = (exp,[((^))+6]) 



p—a 
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It is clear that for a — p one gets g{e) = g{v Ne/a). 
Setting ((A)) = i; we obtain (||). 

Let us summarize the results. Applying the standard 
formalism of maximum entropy principles but formulated 
in terms of nonlinear averages we can explain a large class 
of response functions with asymptotic power law decay, 
including the ubiquitous two-power law. Moreover, using 
the same formalism and standard central limit theorem 
we obtain a universal behavior of fluctuation spectra. In 
both cases the universal laws we propose differ slightly 
from the known ones. To distinguish experimentally be- 
tween them one needs more precise experiments. 

The stay of MC at Universiteit Antwerpen, UIA, was 
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nical University of Clausthal. We are indebted to Prof. 
H. D. Doebner for his comments and his interest in this 
work. MC thanks the Alexander von Humboldt Founda- 
tion for support. 
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